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The gap caused by a strong electric field between the quark surface and nuclear crust of a strange 
star is studied in an improved model including gravity and pressure as well as electrostatic forces. 
The transition from gap to crust is followed in detail. The properties of the gap are investigated for 
a wide range of parameters assuming both color-flavor locked and non color-flavor locked strange 
star cores. The maximally allowed crust density is generally lower than that of neutron drip. Finite 
temperature is shown to increase the gap width, but the effect is significant only at extreme tem¬ 
peratures. Analytical approximations are derived and shown to provide useful fits to the numerical 
results. 
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I. INTRODUCTION 

Strange stars are stars made of absolutely stable quark 
matter Eilli Their existence (i.e., the stability 
of strange quark matter) depends on poorly constrained 
strong interaction properties, and remains to be decided 
by observation or experiment (see @,0 for reviews). If 
strange stars are stable, they contain roughly equal num¬ 
bers of up, down, and strange quarks, but due to the 
higher mass of the s-quark, they will normally contain a 
slight deficit of strange quarks with negative charge, and 
thus have an overall positive quark charge to be compen¬ 
sated by electrons. Even color-flavor locked quark mat¬ 
ter, which is electrically neutral in bulk 0, has an overall 
positive quark charge due to surface effects 111103. The 
quark surface of a strange star is very sharp (the density 
drops from above nuclear matter density to zero within a 
few fm), and since the electrostatic force is weaker than 
the strong force, some of the electrons necessary to cre¬ 
ate an overall charge neutral object will form a thin at¬ 
mosphere with a huge electric field (up to 10^^ V/cm) 
outside the quark phase. This field is capable of sus¬ 
taining a nuclear matter crust decoupled from the quark 
phase by an electron-filled gap. The maximum mass of 
such a crust is approximately 2 x correspond¬ 

ing to the situation where the inner boundary of the crust 
reaches the neutron drip density, 4 x 10^^ g cm“^, where 
neutrons drip out of crust nuclei and get dissolved in the 
quark phase. Smaller crust mass limits would occur if the 
gap is sufficiently narrow to allow direct contact with the 
crust or if the rate of quantum tunneling of nuclei into 
the quark core of the star is large. The crust mass and 
moment of inertia as well as the coupling between core 
and crust play important roles in the understanding of 
strange star properties, and therefore the properties of 
the gap are very important. The width of this gap is the 
main topic of the present investigation. 

The properties of the electron atmosphere and the im¬ 
pact on a nuclear matter crust have been studied by many 

authors, including mis El El El El El El El. These 

studies have involved solutions to the Poisson equation 


for the electrostatic potential with the boundary condi¬ 
tion of electric neutrality deep within the strange star, 
as well as a condition for the potential at large distances. 
For bare strange stars, i.e., stars without a nuclear crust, 
the potential goes to zero at infinity. For strange stars 
accreting nuclear matter on the surface, the condition 
has been taken to be a matching of the electric potential 
from the electrons to the value of the potential required 
in the bulk of the nuclear matter crust (typically several 
MeV, depending on the crust density). This leads to a 
typical gap size (distance from quark surface to the inner 
surface of the nuclear crust) of 10^ — 10^ fm, comparable 
to the distance over which the potential drops by a factor 
of a few. 

The approaches just described do not account however 
for the detailed balance between electrical and gravita¬ 
tional forces and pressure in the transition from gap to 
crust. The potential and its derivative are taken to be 
continuous across the crust boundary while gravitational 
forces and pressure in the crust are assumed insignificant 
(though see El for an approximate inclusion of gravity). 
In this work we expand earlier treatments, investigate 
the effects of including the Newtonian gravitational field, 
and find the detailed structure of the gap and the tran¬ 
sition to the crust. This leads to more narrow gaps than 
previously found and thus constrains the density at the 
base of the crust below the neutron drip density. As the 
temperature increases, the additional thermal reservoir 
of electrons initially widens the gap then narrows it, but 
even for temperatures as high as 10^ K in the gap region 
this effect is small and the width remains constant up to 
temperatures in the range 10^ — 10^^ K depending on the 
crust density. 


II. EQUILIBRIUM OF THE CRUST 

The effective chemical potential of species i, is de¬ 
fined as the change in energy for a unit change in number 
density of species i, while the volume, entropy, and other 
number densities are kept constant. In a steady state as- 
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suming the crust to be isothermal the effective chemical 
potentials of electrons and nuclei should be constant to 
avoid migration. They must therefore be equal to their 
values at the top of the crust, at radius r = R, where 
the electric potential is zero and the chemical potentials 
are equal to the particle masses. Therefore we take the 
effective chemical potentials for electrons and nuclei to 
be 


lj!f =He{r) - e4>e{r) + me4>g{r) =me + nieCpgiR) 

IXe{r) = me + me[(t)g{R) - 4>g{r)] + e4)e{r) 

=I^N{r) + Ze4>e{r) + mN<l>g{r) = rriN + mN4>g{R) 
liN{r) = rriN +mN[(t>g{R) - (t>g{r)] - Ze4>e{r). 

( 1 ) 


For simplicity we describe electrons as well as nuclei as 
Fermi gases characterized by thermodynamic chemical 
potentials fie and /iw- rUe is the electron mass, mw is the 
mass of the dominant nucleus at the base of the crust, 
Z is its charge, is the electric potential and (j)g{r) 

is the gravitational potential for which we will use the 
Newtonian value to get tractable expressions, e^e is typ¬ 
ically of the order of 20 MeV at the strange core surface, 
r = Rs, and (j)g{Rs) = —GMs/Rs ^ —0.2, where Ms is 
the mass of the strange core. 

We note at this point that keeping the effective chem¬ 
ical potentials constant corresponds to hydrostatic equi¬ 
librium which may be seen explicitly by taking the gra¬ 
dient of the above expressions and using the zero tem¬ 
perature identity dPi = nidi^i^ where Pi is the pressure 
and Ui the number density of the i’th component. 
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The total pressure is P = Pe + Pn so the sum then gives 
the usual force balance 


dP d{e(t)e) 
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9_ _ 


= 0 


( 2 ) 


with riq = ZriN — the charge density and p = rriene P 
tuntin the mass density. We show later that Eq. <n is 
convenient when approximating thin crusts. 

To solve for the chemical potentials we take the Lapla- 
cian of Eq. o and use Poisson’s equation in Newtonian 
gravity to get a system of coupled differential equations 
in the chemical potentials and their derivatives. With 
the understanding that the core at r < Rs is comprised 
of quarks and electrons while only electrons are present 
in the gap at Rs < r < Rc, and that nuclei do not ap¬ 
pear until Pn > Sit r > Rc this gives in a compact 


notation 




V = 


dr^ ' Y- dr 
=V^[e(/)e(r) - me(/)g(r)] 

= - 47re^(n+ + Zun - rie) 

— MGrrieipq P TUnTIn P TaeUe) 
d?P n 2 dpN 


( 3 ) 


dr^ j. Pj. 

=V‘^[-Ze4>e{r) -mN(l)g{r)] 
=47re^Z(n^ -1- Ztin — Tie) 

— ATrGrriNiPq + TUnTIn P ^e^e) 


( 4 ) 


Here pq is the approximately constant mass density of 
the quark matter and is its positive charge density 
- we assume that the quark distribution is not signifi¬ 
cantly affected by the crust (see however the recent work 
by Jaikumar et al. 0). In these equations we have used 
rest mass times number density for the matter density in 
Poissons equation. This is a very rough approximation 
for the electrons as they are relativistic except for very 
thin crusts, but in all the calculations below we will ne¬ 
glect the electron contribution to the density completely 
by taking rrie = 0. The densities of electrons and nuclei 
are found from 



fipi) 


1 


( 5 ) 


where p is the particle momentum, pi the statistical 
weight, Ei = the energy and f{p) the Eermi- 

Dirac distribution. Eor a cold Eermi gas with pi = 2 this 
gives 


rii = 


iP - mf)3/2 
Stt^ 


(6) 


while for finite temperatures there will be additional ther¬ 
mal electrons and working in the approximation rrie = 0 
the net electron density (electrons minus positrons) will 
be 


^ “ 37r2 + 3 


( 7 ) 


Equations and 0) can thus be transformed into four 
first order coupled differential equations in the four un¬ 
knowns Pe, Mat, Given appropriate boundary 

conditions at some point - usually either in the charge 
neutral bulk of the strange star core or at the core sur¬ 
face - the system can be solved to find the structure of 
the gap and the transition to the crust. In principle we 
could extend the model to include the entire crust, but 
it should be noted that we have not included the inter¬ 
actions between nuclei (apart from the mean electric po¬ 
tential and gravity) and that the global properties of the 
crust should be described using relativistic gravity and 
a more advanced equation of state - such as the BPS 
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equation of state. Here we restrict ourselves to modeling 
the gap and the lower crust until the point where charge 
neutrality is achieved. The density at this point can then 
be taken as the inner crust density used to parametrize 
the crust models in e.g. [13 


III. APPROXIMATIVE MODELS 

Before we get into the detailed numerical models it is 
worth noticing that a number of illustrative approximate 
models can be calculated analytically, and that as it turns 
out these models reproduce the behavior of the system 
accurately. 

First of all it should be noted that we can get a very 
useful relation from Eq. o by taking the sum 

+ Z^ie{r) = {niN + Zme)[l + (f>g{R) - (f>g{r)] 

( 8 ) 


and its derivative 


d/xjv(r) ^ ^ diie(r) 
dr dr 


= -{rriN + Zrtie) 


d<i>g{r) 

dr 


( 9 ) 


From this and equations (0 0) we can get a qualita¬ 
tive idea of the form of the solutions. For example we 
know that just above the strange core surface the elec¬ 
tric potential should be large and positive and its gra¬ 
dient should be large and negative while the gravita¬ 
tional potential and its gradient can be approximated 
by their Newtonian values. Only electrons are present 
at this point so the charge density is negative, and fie 
will be large and decreasing rapidly but with a positive 
curvature. Since the right hand side of equation 0 is 
almost constant across the gap fiN must increase as jie 
decreases. Thus we also see that jie must decrease mono- 
tonically since jie and jiN would diverge in opposite di¬ 
rections if /ie started increasing while the charge density 
was still negative and /ie(^) was curving upwards. Once 
jiN > m-AT we enter the crust and at some point achieve 
charge neutrality. curves downwards only if the 

charge density is negative and since we expect the crust 
to be charge neutral in bulk must start decreasing 
along with /ig at this point for the crust to remain charge 
neutral. The right hand side of equation 0 is of the or¬ 
der of ^ —10“^^ MeV/fm so for both and to 
be negative they must both be numerically smaller than 
10“^^ MeV/fm. Since will generally be of the order 
of 0.1 MeV/fm close to the surface of the strange star 
core this places high demands on the numerical proce¬ 
dure as it must be accurate across these many orders of 
magnitude. This difficulty stems from the fact that we 
are trying to model the transition from the gap where 
the strong electrical force is dominant to a charge neu¬ 
tral crust held down by the much weaker gravitational 
force - a very stiff set of equations. 

To get something more quantitative than these heuris¬ 
tic arguments we will use the approximation rUg = 0. 


This allows us to solve for /ig in the case of a pure elec¬ 
tron atmosphere with no crust. This can then be used 
to find approximate expressions for the gap width. For 
more realistic solutions it should be taken into account 
that some of the electrons from the crust will spill into 
the gap shielding the strange star core surface charge and 
affecting the gap width. We study the full solutions nu¬ 
merically in a later section, but these are much easier to 
understand (and find) with the insight from the qualita¬ 
tive models below. 


A. The pure electron atmosphere around a bare 
strange star 


For a pure electron atmosphere with rUg = 0 we have 
Me(^) = e0g(r) and in general 

( 10 ) 

Following 0113 we define y = d(e(\)^/dz where 2 : is the 
height above the strange core surface. Assuming a fiat 
strange core surface Eq. 0 then takes the form 


d{e(j)e) 


= 4C-\{e4>er + 7r^T\e4>e)) 


( 11 ) 


where C = y37r/2/e = 5013 MeV fm. This may be 
solved directly for the electric field and potential using 
the fact that both should go to zero at large z: 


- ^ = C-\{e4>ef + (12) 


C 

V27rT 



— sinh ^ 


/ \/27rT y 
(e</>e(-Rs') J 

(13) 


where we give z{e(j)) as this is what we will need later. 
It may be noted that in the limit T ^ 0 equation 
reduces to the expression found by Alcock et al. in 
since sinh~^ x ^ x for x ^ 0 


e(l)e{r) = 


C 


z + C / e(l)e{Rs) 


for T = 0 


(14) 


A plot of this behavior can be found in Fig. ^ Note 
that we have included the case e(j)e{Rs) = 30 MeV. In 
the models below this just allows a crust at neutron drip 
density to remain out of direct contact with the strange 
core surface making it an illustrative case, though the 
choice is difficult to realize for non color-flavor locked 
quark matter equations of state 0 


B. Thin crusts and test particles 

Sufficiently thin crusts should not be expected to 
change the electric potential or the distribution of elec¬ 
trons in the gap from that of a pure electron atmosphere. 
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FIG. 1: Electric potential above the surface of a bare strange 
star. From the top e(j)e{Rs) — 30, 20,10 MeV. We show the 
variation with temperature for e(j)e{Rs) = 10 MeV so from 
the top T = 0,1, 5,10, 20 MeV for these. 


We can then use the solution found in the previous sub¬ 
section and demand that equation m is obeyed at the 
crust boundary. Since we are working in the approxima¬ 
tion rrie = 0 the electronic part of the force balance is 
simply — Ue = 0 which is trivially fulfilled since 

/ie = and dPe = nedjie- For a thin crust we should 
be able to ignore the pressure gradient from the nuclei 
and the remainder of the force balance equation involves 
only the equilibrium between electrostatic repulsion and 
gravitational attraction of positively charged nuclei: 


d{e(j)e) GMsUnttin ^ 

Zun —-^-5-=0 (15) 

dr 

The number density of nuclei then cancels and the equa¬ 
tion simply expresses that each individual nucleus feels 
a balance between the gravitational attraction from the 
strange star core and the electrostatic repulsion from the 
net positive charge below it as described by d{e(\)e)/dr. 
This is the same equation as would be obtained by bal¬ 
ancing a single positively charged test particle above the 
strange star surface. In the case of a cold atmosphere we 
use the derivative of equation for the electric field 
which leads to the relation 

__ GMstun . X 

{Zgap + C/eMRs)f ~ ZRl ' ^ ’ 

or zgap = Rc-Rs- Rs{CZlGMsmNY/^ lO^o fm. 
The gap width of very thin crusts is therefore almost 
macroscopic in size as might be expected since the lower 
layer does not have to support the bulk of the crust. If we 
assume that equation is obeyed identically through¬ 
out the entire crust as it should be if the pressure gra¬ 
dient from the nuclei remains negligible we can solve for 
the potential 


e(t>e{r > Rc) 


GMrriN 


( 17 ) 


where iF is a constant. This potential is a solution to 
Poisson’s equation without a source term, 


V"(e0e) 


(Pjecpe) 

dr‘^ 


^ 2 d{e(l)e) 
r dr 


= 0 


(18) 


showing that the crust is locally charge neutral in this 
approximation. 


C. Cold charge neutral crusts 


In heavy dense crusts the pressure from the nuclei can 
not be ignored, and we would need the precise density 
profile to solve equation The crust would still be 
expected to be charge neutral beyond some point not very 
far from the crust boundary however, so by assuming 
charge neutrality and using equation 0 it is possible 
to estimate fie at the crust boundary. The gap width 
can then be found assuming that the potential in the 
gap can be described by the solution for a pure electron 
atmosphere. 

For a cold crust charge neutrality in the Fermi gas 
approximation is described by the requirement 

37r2 ^ 37r2 ' 

Solving for fie and using relation this gives 


l^e = 


^(miv + Zme)[l +(f>g{R) 




{m% - Z-^Gml){l - 


Z Zrrie) 

{l + cPg{R)-cPg{r))f 



{Z{1 - Z-^G))-\ 

( 20 ) 


Taking rUe = 0 and using (pg (R) — pg (r) 1 this simpli¬ 

fies to 


Me — 


mjv 


(4>g(R) - Mr)) 

-I(Z8/3-l)[0,(i?)-,(,,(r)]2y 


( 21 ) 


Using equation again and p = riNTriN this leads to a 
relation between the size, mass, density and composition 
of the crust 


MR)-Mr) = 


^ fiN — rriN 
rriN 

1/3 


(^8/3 _ 


. 1/2 


= 1.9 X 10“ 


('jLV'Va 


\PdJ 


4/3 


A J 


{z^G _ 1 j 


1/2 


( 22 ) 

(23) 


where pD = 4 x 10^^ g/cm^ is the neutron drip den¬ 
sity. Assuming that charge neutrality occurs shortly af¬ 
ter the crust boundary at Rc^ will not change 
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much between these two points and we can use this 
to estimate jHeiRc)- The crust boundary is defined 
by /iAT = rriN and from equation 0 we then have 
l^e{Rc) = ~ <Pg{Rc)) which Can be used with 

the potential dH to find the gap width 


_ C C 

~^Jie{Rc) eMRs) 


/ n \ / 4 \ 

(w) (s;) 

167 fm 

e0e(7?S')/3OMeV 


=170.5 



- 1/3 

fm - 167 


/ e(l)e{Rs) \ 
V 30 MeV J 



The last equality assumes a crust composed of electrons 
and ^^Fe nuclei. In this approximation the gap width is 
very small near neutron drip density, and for a realistic 
value of e(j)e{Rs) ^ 20 MeV it would be impossible for a 
crust to reach neutron drip density while remaining out 
of contact with the strange star core. If the more strin¬ 
gent limit of Zgap > 200 fm from 0 for the crust to re¬ 
main secure against strong interactions with the strange 
star core is used, the maximum crust density is ^ O.lp^:). 
We will see later that this simple relation describes the 
outcome of the numerical calculations very well. 


D. Charge neutral crust at finite temperatures 

At finite temperatures the electron density is given by 
equation 0, and we will approximate the electric po¬ 
tential in the gap by equation m- As before the crust 
boundary is given by fi]sf{Rc) = tun so iJie{Rc) = = 

^^/S.(j)g by equation 0 with = (f)g{R) - (pgiRc), 

and we will again assume that charge neutrality oc¬ 
curs shortly after the crust boundary so we can solve 
He = ZriN for /S.(j)g. Because of the extra electrons this 
can not be done analytically since the equation for /S.(j)g 
now becomes 

0 =Z(/i^ — — (mAr[l + A(/)^] —^ 

- 3^"(mw[l + A(/)^] - fiN)T^/^Z 


with 


jj^N = {{Sir^p/rriNf^^ + • (26) 

These equations can however be solved numerically for 
A(/)^(p, T) which gives fie{Rc) and Zgap as shown in 
Fig. El Qualitatively we note that changes from the cold 
case are negligible until temperatures around 1 MeV or 
10^^ K are reached - that is until the temperature is com¬ 
parable to pe in the bulk of the crust. For lower densities 
the effect occurs at somewhat lower temperatures, but by 
then the gap is already very large and the temperature 


does not really change this. The bump in the gap size 
seen at high temperatures can be understood from equa¬ 
tions m and (EH). The temperature dependent part of 
equation (1^ shows that and thus e(j)e{Rc) will be 
constant for low temperatures, T pe ^ Rc)i and 
smaller for higher temperatures. Since sinh~^ x ^ x for 
\x\ < 1 and sinh“^ x ^ log2x for x > 1 this can produce 
the bump. Physically the thermal electrons in the crust 
means that pe does not have to be as large to reach a 
certain density, which moves the crust boundary further 
out, while the thermal electrons in the gap screens the 
strange star surface charge and reduces the gap. 


IV. NUMERICAL SOLUTIONS 
A. Numerical procedure and boundary conditions 

To solve equations 0 and 0 numerically we employ a 
shooting method which integrates the equations from the 
quark core to the point of charge neutrality in the crust 
using a standard Runge-Kutta routine. We then vary 
dpeldr at the starting point until equation 0 is fulfilled 
at charge neutrality so the crust remains charge neutral 
beyond this point. Equations 0 would seem to produce 
many more free parameters at the starting point than 
just dpe/dr, but since e(j)e{r) at the quark core depends 
on unknown details of the quark matter equation of state 
we keep it free and explore the solutions for a wide range 
of this parameter. Furthermore it was shown in the pre¬ 
ceding sections that the quantity A0^ = — (j)g{Rc) 

is related to the density of the crust in the charge neutral 
bulk through equations EHli and - and as it turns 
out these relations are fulfilled for the numerical solu¬ 
tions as well. Since dix^/dr may be determined at the 
starting point by equation 0 this leaves dpejdr as the 
only parameter to be determined. With this method we 
can thus choose e(j)e{r) at the starting point (this corre¬ 
sponds to choosing parameters for the equation of state 
for the strange matter core) and a crust - described by 
its density {Acjyg) and temperature at charge neutrality - 
and find the solution to equations @ and 0 describing 
the transition between the core and crust. 

We have so far not specified the exact location of the 
starting point for the integration because the discussion 
above does not depend on this, and because the assump¬ 
tion of a color-flavor locked quark core places the starting 
point at the core surface, whereas for a non color-fiavor 
locked core it should be taken in the charge neutral bulk 
of the core. Color-flavor locked strange quark matter is 
charge neutral in bulk because the pairing minimizes the 
energy for equal quark Fermi momenta leading to equal 
numbers of u-, d-, and s-quarks. However as discussed 
in 1^ when surface effects are taken into account there 
will be a deficit of s-quarks near the surface producing 
an overall positive surface charge. As shown in this 
gives a Coulomb barrier at the surface of the order of 
e(j)e{Rs) ^ 36 MeV. For color-flavor locked stars we thus 
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take the starting point at the surface of the quark core 
and find solutions for different values of e(j)e{Rs) and dif¬ 
ferent crusts. 

In non color-flavor locked stars the deficit of s-quarks 
is global and compensated by electrons in the bulk. The 
surface charge here arises because some electrons near the 
surface leave the core and create an atmosphere outside 
the quark phase. We therefore take the starting point 
“deep” in the core (1000 fm below the surface; our results 
are insensitive to this choice) where we assume charge 
neutrality, n+ ) = ne(/ie, T), 


37r2 37r2 3 


(27) 


The chemical potential for the quark charge de¬ 


fined by this relation, was explored in eg. [J, [111, [1^ 
and found to be of the order of /i+ 25 MeV. However 

the same surface effect which gives rise to the strange 
quark deficit in color-flavor locked stars will reduce the 
number of strange quarks further. We do not explicitly 
include this effect in our calculations. 


B. Chemical composition of the crust 

The chemical composition of the crust depends on its 
origin and the accretion history of the strange star. If the 
crust was created in the supernova along with the strange 
star and has not changed since, it will consist of cold cat¬ 
alyzed matter as described by the BPS equation of state 
in 0- If the crust was accreted from a companion star 
temperatures do not become high enough for the matter 
to reach the equilibrium described by the BPS equation 
of state. Instead hydrogen is burned at the surface to he¬ 
lium which is in turn burned explosively in X-ray bursts 
and leaves a layer of heavy ashes at densities exceed¬ 
ing 10^ g/cm^. The ashes then sink under the weight 
of accreted matter and under the increasing pressure its 
composition changes in a series of electron captures until 
neutron drip sets in. This process has been investigated 
by Haensel & Zdunik, and the chemical composition as 
a function of density was given in [^, assuming that 
the ashes left by the X-ray bursts are composed of ^^Fe. 
The composition of the X-ray burst ashes is a matter of 
some debate however, and as it was shown by Schatz et 
al. in that rp-processes in the X-ray bursts may lead 
to a composition dominated by much heavier nuclei with 
A 106, Haensel & Zdunik recently revised the resulting 
chemical composition in 1^ . However it was also sug¬ 
gested by Schatz et al. that the X-ray burst ashes 
may burn explosively at densities around ^10^ g/cm^ 
powering X-ray superbursts in which case photodisinte¬ 
gration will lead to a composition similar to that found 
in 0. 

For our purpose we chose the relevant nucleus from 
these articles based on the density in the crust at the 
point of charge neutrality - the transition from gap to 


crust takes place over just a few fm so there is no rea¬ 
son to assume any other nucleus will play a role. As 
discussed the actual composition is rather uncertain, so 
we include calculations for each of the three scenarios. 
The charge to mass ratios are very similar, so the only 
significant difference turns out to be related to the dif¬ 
ferent values for the neutron drip density. The chemical 
compositions used are shown at densities below neutron 
drip in Tabled We hereafter refer to the compositions 
in as BPS, HZ1990 and HZ2003 respectively. 

For the HZ1990 and HZ2003 compositions which arise 
from explosive helium burning we assume a composition 
of pure ^He below 10^ g/cm^. This transition should be 
more smooth, but it hardly matters as the gap width is 
already very large at these densities and small corrections 
from the composition will not affect this. 


Ref. [12] 

- BPS 

Ref. [221 - HZ1990 

Ref. [22] - HZ2003 

Pmax 

Nucleus 

Pmax 

Nucleus 

Pmax 

Nucleus 

[g/cm®] 


[g/cm®] 


[g/cm®] 


8.1 X 10® 

56 pg 

1.494 X 10® 

56 Pg 

3.517 X 10® 

“6pd 

2.7 X 10® 

62 Ni 

1.115 X 10“ 

56 Cr 

5.621 X 10® 

“®Ru 

1.2 X 10® 

64 Ni 

7.848 X 10“ 

56 Ti 

2.413 X 10“ 

“6 Mo 

8.2 X 10® 

®4Se 

2.496 X 10“ 

®®Ca 

6.639 X 10“ 

106 Zr 

2.2 X IQio 

®2Ge 

6.110 X 10“ 

56 Ar 

1.455 X 10“ 

106 Sr 

4.8 X 10“ 

®0Zn 



2.774 X 10“ 

106 Kr 

1.6 X 10“ 

78 Ni 



4.811 X 10“ 

“6Se 

1.8 X 10“ 

“Fe 



7.785 X 10“ 

“6Ge 

1.9 X 10“ 

“4 Mo 





2.7 X 10“ 

122 zr 





3.7 X 10“ 

I20sr 





4.3 X 10“ 

“®Kr 






TABLE I: Chemical composition of the crust as a function of 
density for the three scenarios discussed in the text, pmax is 
the maximum density at which a nucleus is present before it 
undergoes electron capture. 


V. RESULTS FOR COLOR-FLAVOR LOCKED 
STRANGE STAR CORES 

We have found solutions for a range of densities, com¬ 
positions and temperatures of the crust as well as for 
different electric potentials at the quark core surface. 
The basic features of one such solution with pcrust = 
4 X 10^^ g/cm^, T = 0, e(j){Rs) = 30 MeV and the crust 
composition in HZ2003 are shown in Figs. mni and El 
The numerical solution conforms with the general be¬ 
havior discussed previously, but it should be noted that 
the actual transition to the crust takes place over less 
than a fm; smaller than the radius of the relevant nu¬ 
clei. A model based on statistical physics should hardly 
be trusted on such scales and the most reasonable con¬ 
clusion would probably be that the density is essentially 
discontinuous at the crust boundary. 

The dependence of the gap width on density, composi¬ 
tion and electric potential at the quark surface is explored 
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FIG. 2: Electric potential and field for a typical solution. 
Note how the field dies out in the charge neutral bulk of the 
crust. The inset focuses on the electric potential. 



FIG. 3: Electron, nucleon and charge density for the solution 
in Fig. 13 The transition to the crust is shown in Fig. 2| 



z [fm] 


FIG. 4: Electron, nucleon and charge density during the tran¬ 
sition to the crust for the solution in Figs. 13 and 13 Nucleons 
appear over a few fm and charge neutrality is attained in the 
crust. 


at zero temperature in Fig. 0 and compared to the ana¬ 
lytical approximation in Eq. (EH). In each case we plot 
(close) to the highest possible density, which is either the 
neutron drip density given in Table ID or the density at 
which Zgap ^ 0. The gap width is insensitive to the choice 
of chemical composition, but it may be noted that only 
for high e(j)e{Rs) is it possible to reach neutron drip den¬ 
sity before the gap width goes to zero. In particular the 
BPS and HZ 1990 compositions are limited by neutron 
drip at e(j)e{Rs) = 30 MeV whereas the gap width goes 
to zero before neutron drip for the HZ2003 composition. 
As noted previously this choice of e(j)e{Rs) is therefore 
quite illustrative and we will often use it below. 



FIG. 5: Gap width dependence on density, crust composi¬ 
tion and electric potential at the quark surface for color-flavor 
locked cores. From left to right e0e(i?s) is 1, 5, 10, 20 and 
30 MeV. Gurves for different composition are almost indis¬ 
tinguishable. The analytical approximations are based on 
Eq. (1^ using the crust composition from the BPS equation 
of state. The dashed line shows the solution of Eq. (1^ with 
Zgap = 0 - i.e. the density at which the gap width goes to 
zero for a specific e(j)e{Rs)- 

Figure 13 shows the temperature dependence of Zgap 
with e(j)e{Rs) = 30 MeV. Temperatures up to 100 MeV 
are shown to illustrate the full range of solutions, but one 
should keep in mind that the Coulomb barrier is only 30 
MeV or less, and so can not hold the nuclei at these tem¬ 
peratures. e(j)e{Rs) is set very high to allow high den¬ 
sities, and for simplicity we only vary the composition 
for the maximum density curves discussed above. We 
recover the qualitative features of the analytical approx¬ 
imation. Except perhaps for the very early stages after 
formation in a supernova explosion crust temperatures 
of isolated strange stars are expected to be much smaller 
than 0.1 — 1 MeV where the temperature effects for the 
gap width become noticeable. Interesting temperatures 
may be reached in accreting binary systems to which we 
return in a later section. It may also be relevant that 
the increase in gap width allows a strange star to sus¬ 
tain a crust shortly after its formation in the supernova 
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even when temperatures are very high - whether a crust 
would actually form at that stage is a different matter. 


10000E 


1000 ,; 


100 , 


— Analytical approximation 

o-€.BPS,p=Pn,ip 


0-0 BPS, P=pDrip''l0 

^BPS,p=PQ,,p/100 

<^BPS,p=p^.A000 

V—S' HZ2003, p=6.5 ♦ lo" g/cm^ 
>_^HZ1990.p=p^,p 

A A A -A 




J—e—0 ® > 



0.001 


0.1 1 
T [MeV] 


100 


FIG. 6: Gap width dependence on temperature, density and 
crust composition for color-flavor locked cores with e(l)e{Rs) = 
30 MeV. Densities increase from top to bottom and curves 
without points show analytical approximations based on 
Eq. (1251 . For simplicity we only vary the composition for 
the maximum density crusts. The HZ2003 curve has a very 
low gap width because its density is not limited by neutron 
drip. 

Perhaps a better measure of the stability of the crust 
is the transmission coefficient for ions through the gap, 
r, which following Alcock et al. 0, can be found in the 
WKB approximation by 


le+06 


10000 - 


100 


0.01 



.1-r-nrnrrTTTTi-r-rnmrrTT 

1-1—r T-T-TTTr|-1—n-i-rmi-i i i 


- 



- 


<3— ^^^^^ 



1 

0—0 BPS, -log x(p) 



- 

M HZ2003, -log T(p) 




^HZ1990, -logx(p) 




<3—< BPS, -logT(T) 

> [> 



v-v HZ2003, -log x(T) 




^ HZ1990, -log t(T) 



- 

- - Break even 


- 


— Analytical approximation 



1 

1 ■ . ■M..I■ ■ . ■■■■■1■ ■ . MM. 

[ 

1■ . ■ M.l.l ■ . ■ ...mI.. 

] 1 


le+06 


le+07 


le+08 


le+09 


le+10 


T [K] or p [g/cm ] 


le+11 


le+12 


FIG. 7: Dependence of the transmission coefficient on den¬ 
sity at T = 0, and on temperature at maximum density with 
e(j)e{Rs) — 30 MeV. The kink in the HZ2003 density depen¬ 
dence is at the assumed transition to a hydrogen crust. The 
horizontal line marks the break even estimate for accreting 
stars in the text. 


bering that e(j)e{Rc) = this gives 


log(r) = 

2V2CZ _i / mivA(/)g 

\ZeURs)) 


_ mN^(f>g 
Ze4)e{Rs ) 
Ze(f)e{Rs) 

mN^4>g 


1/2- 


(32) 


r = exp 


pZ = Zg 

—2 / \k\dz 

Jz=0 


(28) 


where k = is the wave number, and the 

chemical potentials are known from the numerical solu¬ 
tions. A few examples are shown in Fig. [3 Here we again 
choose e(j)e{Rs) = 30 MeV and plot the dependence on 
density at zero temperature, and the dependence on tem¬ 
perature for the maximum density crust of each compo¬ 
sition. 

The integral in Eq. (123 can be approximated at T = 0 
by 


logr ~ 


rzg / 2 ; 

-2V2mi\j / (-e())e(r) — A(pg ) dz (29) 

Jo \^N J 

px{Rc) 

-2V2CZ - A(j)gY/^dx, (30) 

J x(Rf^) 


where 




X = I - e4>e{r) 

\mN 


1/2 


(31) 


using Eq. d) for fiN and Eq. HI for e(j)e{r). Remem- 


where can be approximated by Eq. (I 23 . As seen 
from Fig. □ Eq. (EH) is a very reasonable approximation. 
Assuming the HZ2003 crust composition at drip density, 
the low density limit of Eq. (EH) is given as 


log(r) 


-2.5 X 10'‘ 


TT 

2 



1.84 


(33) 


The transmission coefficient may be better appreciated 
by estimating its value when transmission of ions into 
the core through the barrier breaks even with typical 
accretion rates in binary systems. That is, we require 
the balance: 


- < Nion X f XT , (34) 

tun 

where M is the mass accretion rate onto the strange star 
from a companion, Nion is the number of ions, whose mo¬ 
tion about their lattice position allows them to strike the 
barrier, and / is the oscillation frequency of this motion. 
Nion iiiay be taken as roughly the number of ions within 
one lattice distance, a ^ 200 fm, from the crust bound¬ 
ary, and following again Alcock et al. 0, the oscillation 
frequency should be less than 1 MeV. For a strange star 







































9 


accreting 10 ^^MQyr ^ we then arrive at the condition 

- logr > 44.6 - log + lo* 

( 4.3xl0ng/cm» ) +‘“® 

(35) 

as indicated in Fig. □ (all logarithms are base e). It may 
then be noted that the maximum density models shown 
in Fig. [7|- except the one using the BPS chemical com¬ 
position - will interact with the core at low temperature 
and are not stable at these densities. Only the BPS equa¬ 
tion of state is thus able to remain stable at neutron drip 
density even with e(j)e{Rs) = 30 MeV. For other compo¬ 
sitions and lower e(j)e{Rs) stable crust densities are lower 
depending on the temperature, see Fig. im 



FIG. 9: Electron, nucleon, quark charge density and total 
charge density for the same solution as in Fig.|3 The transi¬ 
tion to the crust is shown in Fig. cni 


VI. RESULTS FOR NON COLOR-FLAVOR 
LOCKED STRANGE STAR CORES 

We have performed similar calculations for the case of 
non color-flavor locked cores (“ordinary” strange stars), 
and the structure of a solution with pcrust = 4 x 
10^^ g/cm^, T = 0, /i+ = 30 MeV and the HZ2003 crust 
composition is shown in Figs.|Hlini andimi Charge neu¬ 
trality is imposed at the starting point 1000 fm into the 
core and as we approach the surface electrons become 
scarcer and consequently the charge density and electric 
field rises. At the surface the quarks and their posi¬ 
tive charge vanish, the charge density becomes negative 
and the field starts decreasing. The solution from the 
surface out looks like the color flavor-locked case with 
e(/)e{Rs) — pe{Rs)- 



-1000 -800 -600 -400 -200 0 


z [fm] 

FIG. 8: Electric potential and field for a solution in the non 
color-flavor locked case. The inset focuses on the electrical 
potential. 
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EIG. 10: Electron, nucleon, quark charge and total charge 
density during the transition from gap to crust for the same 
solution as in Eigs. |H|and|ni 


We again calculate the dependence of the gap width on 
density shown in Fig. T = 0 and different choices 

of = 10, 20, 30, 40 MeV from left to right. We 
only show models using the HZ2003 crust composition, 
since the gap width is in this case so insensitive to the 
choice of composition that the curves would be indistin¬ 
guishable. Note again that only for the highest choice of 
= 40 MeV does the crust reach neutron drip density 
before the gap goes to zero. The gap width is approxi¬ 
mated by Eq. (1^ using e(j)e{Rs) = = 10 MeV. For 

a bare strange star we would have e(j)e{Rs) = 3/i+ /4 Q, 
whereas for dense crusts e(j)e{Rs) — ^ but since large 

gap widths are insensitive to e(j)e{Rs) we use e(j)e{Rs) = 
/i+ in all analytical approximations. 

The temperature dependence of the gap width is ex¬ 
plored in Fig. El for the case of = 30 MeV with 
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FIG. 11: Dependence of the gap width on density using 
the HZ2003 crust composition. From left to right /xj = 
10, 20, 30, 40 MeV. The analytical approximation explained 
in the text assumes /xj = 10 MeV. 


the HZ2003 crust composition. We note that the gap 
width declines steeply at T ^ 10 MeV and reaches zero 
at T ^30 MeV in contrast to the color-flavor locked 
case. This behavior is caused by the condition of charge 
neutrality at the starting point in the core (Eq. (I27|) b 
which reduces /Xe significantly at temperatures compara¬ 
ble to /x+. Other that that we note again the bump in the 
gap width caused by the condition of charge neutrality 
in the crust. The lowering of /Xe entirely suppresses the 
bump at high densities. The gap width is approximated 
in Fig.Elby Eq. (tT!^ using e(j)e{Rs) = = —1000 fm) 

found from Eq. e(j)e{Rc) is estimated as before from 
Eq. 0 and 



FIG. 12: Gap width dependence on temperature and density 
for non color-flavor locked cores with /xj = 30 MeV. The 
maximum density is taken to be the density at which Zgap — 1 
fm as this is as close to zero gap width as is numerically 
reasonable - note that variations in the numerical results are 
significant only at this scale. The analytical estimate assumes 
P — Pmax /10. 


The variation of transmission coefficient with density 
and temperature is shown in Fig. ^1 We again note 
that the bump is suppressed compared to the color-flavor 
locked case and that only crusts with the BPS composi¬ 
tion can reach neutron drip before they become unstable 
against transmission through the Coulomb barrier. 



FIG. 13: Dependence of the transmission coefficient on den¬ 
sity at zero temperature and temperature at maximum den¬ 
sity for non color-flavor locked strange star cores with /xj = 
30 MeV. The kink in the HZ 1990 and HZ2003 dependence 
corresponds to the assumed transition to a hydrogen crust. 
The break even estimate is explained in the text and the an¬ 
alytical approximation is based on Eq. (1321 with the HZ2003 
crust composition at drip density and e0e(i?s') = Pq- 


VII. DISCUSSION 

We have expanded the treatment of gaps below strange 
star crusts to include effects of pressure and gravity on 
the crust nuclei and found the structure of the transition 
from gap to crust. Overall our results are as one would 
expect with a very sharp transition from gap to crust, gap 
widths ranging from a few fm at densities near neutron 
drip over a few thousand fm at densities close to the 
maximum to 10^^ fm for very thin crusts. Crust densities 
are consequently limited below neutron drip unless one 
assumes an extreme Coulomb barrier height. If the gap 
width must be at least a crust lattice distance 200 
fm) and a reasonable value of e(j)e{Rs) — 20MeV (or 
e(j)e{z = —1000 fm) ^ 20 MeV) is assumed then the crust 
density is limited to a few times 10^^ g/cm^. 

The variation of the gap width with temperature is 
noteworthy in that it initially increases with tempera¬ 
ture which is a qualitatively new feature. However this 
increase only takes place at temperatures below 1 MeV 
(10^^ K) for crust densities which give very small gap 
widths of the order of 1 fm such as in the case of the 
HZ2003 crust at maximum density in Eigs. IHI and 0 A 
realistic crust should have a sufficient gap width to be 
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stable against strong interaction with the core, and exam¬ 
ples of such crusts could be the HZ 1990 and BPS crusts 
in Figs. 0 and Q - the HZ1990 crust is almost stable, 
and lowering the density a little to make it so would not 
change the dependence on temperature much. If we ap¬ 
ply the criterion in Eq. for crust stability we can 
find the highest stable density at given temperature and 
Coulomb barrier height, e(j)e{Rs)- Numerical results for 
this are shown in Fig. IT^and we see that temperature ef¬ 
fects are then not important below 10^^ K. If we demand 
at least a crust lattice distance in non accreting systems 
the crust becomes even more stable against variations in 
the temperature. 



T[K] 


FIG. 14: Highest stable density according to Eq. dD at given 
temperature and e(l)e{Rs) = 1,5,10,20,30 MeV from below. 
The crust composition is from HZ2003 and we stop at the 
temperature corresponding to neutron drip density. 

The thermal structure of accreting strange stars was 
investigated by Miralda-Escude et al. who found 
temperatures around 10^ K at the crust boundary for 
strange stars accreting 10“^^ to 10 “^-^^ Mq yr“^ assum¬ 
ing that 10 MeV is released per nucleon converted to 
strange quark matter (bag constant H = 103 MeV/fm^). 
In this case the gap width would not be affected at all, 
and one might as well use the zero temperature results. 
Superbursters however are known to have much higher 
accretion rates around 0.1-0.3 times the Eddington ac¬ 
cretion rate, M^dd ^ 10 “^ Mq yr“^, and in a recent in¬ 
vestigation Page and Gumming found temperatures 
for such stars in the range 10^ —10^ K at the crust bound¬ 
ary. This work was motivated by the difficulties in achiev¬ 
ing carbon ignition at observed column depths in super¬ 
burst models on neutron stars - see Gumming et al. 2005 
for a recent review - and it was found that super¬ 
bursts may ignite at the right column density on strange 
stars provided that neutrino emission in the core is slower 


than direct Urea which implies that strange quark matter 
should be a color superconductor. During the superburst 
itself temperatures in the range 2 — 7 x 10^ K are reached 
for a few hours This is again too low for signif¬ 

icant temperature effects and for stable crusts one may 
as well use the zero temperature expressions. 

The present treatment could be improved by includ¬ 
ing general relativity instead of Newtonian gravity, and 
by using a less simplistic model for the equation of state 
for the crust. Such improvements are necessary to calcu¬ 
late the properties of the whole star, but we expect the 
changes in the gap properties to be minor compared to 
the results presented here. The same can be said about 
a more detailed treatment of the quark phase, which has 
here merely been treated as a source of a given total mass 
and starting value for the electron chemical potential. 
However, as pointed out in , the electric potential in 
the outer layers of the strange core could be qualitatively 
different from the behavior assumed here if a mixed phase 
of quark nuggets and electrons can be realized. 

Another important oversimplification involves the as¬ 
sumption of a constant pairing gap all the way to the 
quark core surface in the treatment of quark matter in 
the GFL-phase. In fact the pairing gap Acfl may be 
envisaged to vanish within a surface distance of order 
1 /Acfl ~ 2-20 fm for gap energies between 100 and 
10 MeV. Possibly a range of phases appear in the extreme 
uppermost layers, including the possibility that the very 
surface resembles the ungapped ordinary quark phase 
also discussed above. A detailed treatment of gapped 
phases near a surface remains to be performed. 

In conclusion the results of our expanded treatment, 
while quantitatively somewhat different, are qualitatively 
consistent with previous work in the literature with re¬ 
gards to gap width, transmission coefficient and possible 
crust densities. The increase in gap width with tempera¬ 
ture is new but significant only at very high temperatures 
not likely realized. Perhaps the most useful results of the 
investigation are the various analytical approximations 
which have been derived and shown to fit the full numer¬ 
ical solutions very well, since these provide better phys¬ 
ical insight and can be easily used in models for strange 
star rotation, glitches, accretion and instabilities. 
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